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ABSTRACT 

The Clouds and Earth Radiant Energy System (CERES) scanning radiometer was designed to provide high 
accuracy measurements of the radiances from the earth. Calibration testing of the instruments showed the presence 
of an undesired slow transient in the measurements of all channels at 1% to 2% of the signal. Analysis of the 
data showed that the transient consists of a single linear mode. The characteristic time ul this mode is 0.3 to 
0.4 s and is much greater than that the 8- 10-ms response lime of the detector, so that it is well separated from 
the detector response. A numerical filter was designed for the removal ot this transient from the measurements. 
Results show no trace remaining of the transient after application of the numerical filler.. The characterization 
of the slow mode on ihe basis of ground calibration data is discussed and flight results are shown for the CERES 
instruments aboard the Tropical Rainfall Measurement Mission and Terra spacecraft. The primary influence ot 
the slow mode is in the calibration of the instrument and the in-flight validation of the calibration. This method 
may be applicable to other radiometers that are striving for high accuracy and encounter a slow spurious mode, 
regardless of the underlying physics. 


L Introduction 

The Clouds and Earth Radiant Energy System (CE- 
RES) scanning radiometer (Barkstrom 1990; Wielicki 
et ai. 1997) was designed to prov ide high accuracy mea- 
surements of the radiances from the earth in three wave- 
length bands: total (0.2-100 /cm), shortwave (0. 2-5.0 
/cm), and the longwave window channel (8-12 /an). 
From these broadband radiance measurements, the ra- 
diative fluxes leaving the "top of the atmosphere" can 
be computed. CERES instruments flew aboard the Trop- 
ical Rainfall Measurement Mission (TRMM) and the 
Terra spacecraft and are to fly aboard the Aqua space- 
craft. The TRMM was placed into orbit in November 
1997. Its orbit covers 35°N-35 U S and precesses through 
all local times every 48 days. In addition to a CERES 
instrument, the TRMM spacecraft includes a precipi- 
tation radar, a microwave imager, a Visible and Infrared 
Scanner and a Lightning Sensor. The PR permits eom- 
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putation of latent energy release within clouds as a con- 
comitant to the radiant flux at the top of the atmosphere 
in order tt> help understand the energetics of the Tropics. 
The Terra spacecraft was launched into a sun-synchro- 
nous orbit crossing the equator at 1030 LT in December 
1999. Terra includes the Moderate Resolution Imaging 
Spectroradiometer (MODIS). which provides cloud data 
to help diagnose cloud properties to aid in understanding 
the CERES radiances. The Aqua spacecraft will also 
carry a MODIS and infrared and microwave sounders 
for retrieving atmospheric temperature and humidity 
profiles. Its orbit will be sun synchronous, crossing the 
equator in midafternoon. 

The CERES instruments are calibrated in vacuum to 
better than 1 c /v accuracy and precision approaching 0.3% 
(Lee et al. 1996, 1998, 1999, 2000; Barkstrom et al. 
2000; Smith et al. 2000). As a consequence, it is possible 
to detect small imperfections in the dynamic responses 
of these instruments. Calibration testing ot the proto- 
flight model (PFM) of the instrument, which flew on 
the TRMM spacecraft, showed the presence of an un- 
desired slow transient in the measurements from all 
three channels. This transient is about l%-29r of the 
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SAMPLE NUMBER 

Fk;. F Response of the total channel of the PFM as function ol 
sample number as it views internal hlackbody. 

true signal, depending on the channel, and damps out 
over a characteristic time of 0.3-0. 4 s. At the time the 
transient was found, it would have been prohibitively 
expensive to modify the hardware to eliminate the tran- 
sient, and there is no guarantee that such an effort would 
have succeeded. 

One of the major effects of the transient is that it 
corrupts the measurement when the instrument scans to 
stare at space for a zero radiance value (the “space 
look"), causing an offset in all measurements during 
the preceding and succeeding scans. Another major ef- 
fect is that during in-flight calibration the measurement 
does not approach its asymptotic value during the time 
available, thereby creating errors in the computed gain 
of the channel if not taken into account. Also, as the 
instrument scans from the space look across the earth, 
the transient causes the measurement to be low as it 
scans onto the limb. Finally, as the instrument scans 
from a cold bright cloud to a warm dark surface, the 
transient w ill cause a loss in sharpness of the change. 

In this paper we discuss a numerical filter for the 
removal of this transient from the measurements. In or- 
der to remove this spurious signal, it is necessary to 
characterize its behavior, which is defined by two con- 
stants. A method for deriving these constants from the 
data is developed. Finally, results of application of this 
numerical filter to the CERES instrument aboard the 
TRMM and Terra spacecraft data are shown. No trace 
of the transient is discernible after application of the 
numerical filter. For the in-flight calibration results quot- 



Fk;. 2. Asymptotic portion of response of total channel ot PFM as 
function of sample number as it views internal hlackbody. 


ed by Lee et al. (1998, 2000), the numerical filter was 
used. 

The slow response was found because the calibration 
permitted accuracies of a fraction of a percent. At this 
level, one expects many problems to occur. The present 
technique may be applicable to other instruments which 
encounter response transients in operation. 

2. Background 

The CERES instruments each have three channels: a 
total channel (0.2-50+ jum), a shortwave channel (0.2- 
5.0 /Am), and a window channel (8-12 gm). These in- 
struments w ere calibrated in the Radiometric Calibration 
Facility (Lee et al. 1996). Figure 1 shows a time history 
of the output of the total channel of the PFM as the 
instrument scans onto and stares at the wide field-of- 
view hlackbody in the Radiometric Calibration Facility. 
The output rises very quickly to about 98% ot its peak 
value, then has a slow transient in its approach to its 
asymptotic value, shown in Fig. 2. The design response 
time of the instrument — that is, the primary thermal 
mode of the detector — is on the order of 10 ms, which 
is the data-sampling rate. This mode is shown by the 
initial rapid rise. The slow transient that follows can be 
characterized as a slow mode with a characteristic time 
on the order of 300— 400 ms and an amplitude of 1 %- 
2 % of the signal. Similar results were found for the 
longwave window channel (8-12 gm) and the short- 
wave channel (0.2-5. 0 jum). In order to treat the slow' 
mode, it is first useful to examine the mechanism of the 
slow mode and to consider the various time scales of 
importance to the CERES instrument. 

a. Mechanism of slow mode 

Each channel of the CERES instrument consists of a 
Cassegrain telescope, which focuses radiation from the 
scene onto a thermistor bolometer (Fig. 3). This active 
thermistor bolometer is mounted on a heat-sink disk. 
The active thermistor bolometer is connected with a 
bridge circuit to a compensator thermistor bolometer, 
which is mounted on the back side of a second heat- 
sink disk in an enclosure in order to maintain constant 
radiation and temperature. The active and compensator 
thermistor bolometers are selected from those available 
to be as nearly identical as possible. The two disks are 
bolted together w ith an indium gasket between them to 
increase thermal conductivity. A titanium gasket reduces 
thermal conduction from the heat-sink disks to the field 
stop. In essence, the measurement is proportional to the 
difference of temperature between the active and com- 
pensator thermistor bolometers. The temperature of the 
disk subassembly is maintained as nearly as possible at 
a constant level by an electric feedback loop, which has 
a time response of several minutes. 

The slow mode has been simulated quite closely by 
a highly detailed numerical model (Priestley 1997: 
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Fici. 3. Design of radiometer. 


Haeffelin et al. 1997). The spurious mode is due to 
heating of the heat-sink disk to which the active therm- 
istor bolometer is attached, as indicated in Fig. 4. When 
radiation from the scene is focused onto the active 
thermistor bolometer, it heats up with a characteristic 
time of 8-10 ms to provide the measurement. As the 
radiation continues, the heat-sink disk to which the ac- 
tive thermistor bolometer is attached heats up, permit- 
ting the active thermistor bolometer temperature to rise 
and increase the temperature difference between it and 
the compensator thermistor bolometer and increasing 
the measurement. The characteristic time for this mode 
of heating is 0.3-0.4 s. 

b . Timescales 

The timescales relevant to the CERES scanning ra- 
diometer are showm in Fig. 5. The output is sampled 
every 10 ms. Any time response of the radiometer faster 
than 10 ms cannot be resolved from the data and must 
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Fici. 5 . Characteristic times of radiometer. 


be included in the point response function. The radi- 
ometer looks at space every 3 s in order to get a zero 
radiance reading. The zero radiance reading is linearly 
interpolated between space looks. Over a period of sev- 
eral minutes the zero radiance value may drift signifi- 
cantly. Also, the temperature-stabilization feedback 
loop for the disk subassembly acts over a period of 
several minutes. Any time response longer than the time 
between space looks will be attenuated by the subtrac- 
tion of the zero radiance readings from the measure- 
ments. Between the data-sample period and the space 
look period is a 2 1/2 order of magnitude span within 
which one hopes there are no detectable time responses. 
However, in the CERES scanning radiometer this spu- 
rious mode appears with a 0.3-0.4-s time response at 
the 2% level. This mode is sufficiently slow that its 
effects can be uncoupled from the point response func- 
tion and resolved by the data, so that it is possible to 
remove any significant effects by a numerical filter. 

3. Numerical filter formulation 

a. Formulation of slow mode response 

Figure 6 shows the basic concept upon which the 
numerical filter is based. As the result of an input ra- 
diance r(t), a signal u(t) is generated. The signal can be 
expressed in terms of an integral over the field-of-view 
of the scene radiance weighted by the point response 
function. The angular distance in the scan direction is 


Scene 

Energy 



Fig. 4. Mechanism causing slow mode. 
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h(dldt) are linear operators consisting of polynomials in 
the time derivative. The degree of the g polynomial is 
the number of significant modes which contribute to the 
transient response and the degree of h is one less. Equa- 
tions (2) and (3) can be combined to give the spurious 
transient due to the slow modes in terms of the mea- 
surement: 


Fig. 6. Block diagram of radiometer response. 

mapped into time t and the variation of the scene in the 
direction normal to the scan is ignored in the present 
treatment, reducing the problem to one dimension. Thus, 
the signal is 

u(t) - I P(t' ~ t)r(t') dt\ ( I ) 

J K)V 

The data rate, instrument field-of-view (FOV) and the 
signal conditioner (a four-pole Bessel filter) have been 
matched to the time response of the detector so as to 
minimize alias and blur errors in the retrieved field 
(Manalo and Smith 1991; Manalo et al. 1991). 

A basic principle of design is that the signal formation 
should not pass significant power at frequencies higher 
than the Nyquist frequency. All frequency components 
in the scene r(t) faster than the Nyquist frequency are 
attenuated from u(t) by the sensor design, which in- 
cludes selection of the field-of-view size and shape, the 
detector response time and the Bessel filter. These high 
frequency components are a major distinction between 
u(t) and r(/). The second major distinction is that the 
spectral components of u(t) with frequencies near the 
Nyquist frequency are reduced in amplitude, which is 
blur, and usually delayed in phase, which is distortion. 
The phase delay is accounted for by locating the data 
at the centroid of the point response function. 

Unfortunately, due to heating of the heat sink of the 
active thermistor-bolometer, one or more spurious 
modes are superimposed on the signal. The response of 
the spurious modes as they appear in the output is de- 
noted as u(f). The measurement w(t) is thus given by 
the detector output with the spurious modes superim- 
posed: 

vv(r) = u(t) + i At). (2) 

In order to achieve the desired accuracy from the mea- 
surements vv(f), it is necessary to remove the spurious 
mode v(t). We wish to retrieve the undistorted signal 
u(t\ for which the design is already optimized on the 
basis of alias and blur errors, and not r(t). 

The spurious modes are assumed to be driven by heat 
going into the sensor, which is proportional to the un- 
distorted signal «(/), so that (see appendix) 



Although the present problem is due to transient heat 
flux between the active and compensating thermistor 
bolometers, Eq. (3) does not depend on the details of 
the underlying physics, but only states that the process 
is governed by a linear differential equation. Also, the 
equation has constant coefficients, w'hich simplifies the 
present application so that the numerical filter likewise 
has constant coefficients. A system with time-varying 
coefficients could be treated in a similar manner. 

At this point a number of modes could be included; 
the a priori assumption is made that only one mode is 
needed to describe the spurious transient of the CERES 
scanning radiometer. The justification (or rejection) of 
this assumption wall depend on the results. Thus the 
operator on the left hand side of Eq. (3) reduces to first 
order: 



where A is the inverse of the characteristic time of the 
mode. The operator on the right-hand side is then of 
0th order and reduces to cA, where c is the loading of 
the slow mode into the signal. For a single mode spu- 
rious response, Eq. (4) reduces to 

dv ... 

— + A(1 f c)v — cAvr. (6) 

dt 

The solution of Eq. (6) can be expressed in terms of 
the measurements w(r) as 

i ;(/) = c\ J Hi - ?')»•(/'} lit'. (7) 

where /(/) is the impulse response function for the op- 
erator (g + c*A). In particular, for the single-mode case, 

v{() = cX j expl — A( 1 + c)(t - dt\ (8) 

Equation (8) gives the effect of the spurious mode on 
the measurement in terms of the measurement history. 


g 



v{t) = h — ]w(f). 
\dtj 


(3) 


Because conduction is a linear process, gidldt) and 


b. Derivation of numerical filter 

The measurement u(T) is known only at discrete 
points Because the spurious mode is slow, vv(T) can 
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be approximated as piecewise constants based on vv(/ ; ). 
Equation (8) may thus be written as 

v(k) = /;, ph Ar (/), (9) 

i ■ i 

where 

p {l ~ expj - AA/(1 + c)l 

p ] = c( 1 - /?„)/( I + cl ( 10) 

Equation (9) is a very inefficient form for application. 
This solution can be generated recursively by writing 
Eq. (9) for v(k — 1 ), whence it is seen that 

v(k) = p t) v(k ~ 1) + p t w(k). (11) 

The tilde is used to denote that this is a computed es- 
timate of the quantity. Equation (11) can also be derived 
as the solution of the difference equation which cor- 
responds to the differential equation (6) in finite form. 
Equation ( 1 1 ) permits the efficient computation of the 
spurious mode recursively during data processing. Once 
the 77(A) is computed, the retrieved measurement can be 
computed as 

Ti{k) = w(k) - 77(A). (12) 


c. Applic ation of numerical filter 

The modal loading factor c and the modal time con- 
stant A 1 are different for each detector due to manu- 
facturing differences. Once these constants are deter- 
mined from calibration data, the /^ £> and /?, are computed. 
For each measurement ve(A), the spurious mode ampli- 
tude 77(A) is computed using Eq. (1 1) and the retrieved 
measurement 77(A) is computed by use of Eq. (12). The 
algorithm requires only the storage of the immediately 
preceding transient amplitude v(k — 1). 

In order to start the filler at the very beginning or 
after a disruption of data How, the transient mode is 
assumed to be in equilibrium with vv(A) for expediency, 
so that we let 


for the initial conditions only. Equations (II) and (12) 
are then used for processing the data. After a few char- 
acteristic times A 1 , the effect of the initial condition 
will be damped out exponentially. Standard procedure 
is to begin count conversion to raw radiances with the 
scanning of space. During this time, the initial condition 
will damp out. 

d. Response of measurement to step input 

The response of the measurement to a step radiance 
input is now considered. The response of the undistorted 
output u(t) as the radiometer scans onto a unit step input 


is denoted as Fit), where Fit) — » 1 for large t. The 
spurious mode is then given from the solution of Eqs. 
(3) and (4) as 

v(t) - Gc\ j exp|-A(f - f))F(t r ) dt‘ . (13) 

Even though F(t) is a smooth function at the timescale 
of the pixel response during calibration, at the timescale 
of the slow mode it can be considered to be a step 
function. In order to reduce errors due to this approx- 
imation, the step is assumed to occur at the centroid of 


the pixel, which we denote as For t > /,, F{t) ~ 1 
and v(t) can be written as 

v(t) = GY‘| 1 — e Mt (14) 

The measurement is then 

wit) = G -F Cell — e Mt ' !> f (15) 

and the asymptotic value of the measurement is 

= G(1 + c). (16) 


Equation (16) show s that the slow' mode does not vanish 
at large /, but takes a steady-state value. 

The slow mode affects the gain computation of the 
instrument. In the absence of the slow mode, the in- 
strument counts for a unit input would be q — G. During 
ground calibration, the total and longwave window 
channels stare at the blackbody for more than 3 min. 
During this time, the slow mode does not vanish, but 
approaches a nonzero steady state Gc. Equation (16) 
shows that the instrument counts are q — G(1 T c). The 
apparent gain G' is computed on the basis of this steady 
state, thus is G' = ql( 1 + c) and includes the steady- 
state value of the slow mode. This apparent gain G' is 
reduced from the true gain G: 

G f = G/( 1 + c). (17) 

In order to use the apparent gain, that is, the result 
from the steady state, the factor (1 + <■) is included in 
Eq. (12): 

77(A) = rvv(A) - 77(A)|(t + r). (18) 

For the CERES computations, Eqs. (II) and ( 1 8) should 
be applied to filtered radiances, after accounting for the 
offsets and space clamp. This is most easily done by 
incorporating the factor (I + c) in the p {] and p ] . By 
w riting the equations in this manner, the numerical filter 
results match the results obtained from the steady-state 
calibrations. One advantage of this procedure is that the 
steady state results using the numerical filter are inde- 
pendent of the modal amplitude c. In the data process- 
ing, after the instrument counts are converted to radi- 
ances, the numerical filler is used to retrieve the undis- 
torted radiance. 
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I k;. 7. ln|(n lsv \\V))fiu ~ u m .J| ns function of time as total 
channel scans to view internal blackboily. 


4. Determination of slow mode characteristics 
from ground calibration data 

In order to apply the numerical filler to eliminate the 
transient mode from the data, it is necessary to deter- 
mine the two parameters A and <\ which govern the 
dynamic response of the slow mode, from the ground 
calibration data. In ground calibration testing the slow 
mode is most evident when the total or longwave win- 
dow' channel scans from the cold reference (Lee et al. 
19%) to stare at the wide-tield-of-view (WFOV) black- 
body for a period of more than 3 s. Similar data are 
also taken when these channels stare at the internal 
blackbodies for the purpose of calibrating the internal 
blackbodies. Scanning of the total channel from the cold 
reference to stare at the WFOV blackbody or the internal 
blackbody approximates a step input to the sensor. From 
Eq. (16) the step response of the system can be ex- 
pressed in terms of counts as 


wit) = vr ds . 


(vr a 



(19) 


From this form we can write 


vr asv 

- w(t) 





- ln(r) 


w\„ v 

- w rnjn 



in( 1 

+ cl 


(20) 


A plot of ln((u\ ls> . - u*( / ) )/( w. 1S j. " H’ lllln )| wall thus have 
a slope of -A when the basic time response of the 
detector has damped out. leaving only the slow' mode 
response. Figure 7 show's ln|(u a ^ - vr(/))/fn , dS> . — vv ini J] 
as a function of time for the total channel of the CERES 
PFM as it view's the internal blackbody during ground 
calibration testing. The plot shows three regions of re- 
sponse. In region I the instrument is scanning onto the 
calibration source and the detector and signal condi- 
tioning circuitry are coming into equilibrium. Because 
the input is not a perfect step function the instrument 
response is quite complicated in this region. In region 
II the response is linear, indicating that the only variation 
is due to the slow' mode and Eq. (20) applies. The slope 
of the response in this region is - A. the inverse of the 
characteristic time. In region III the slow' mode is small 


enough that the noise begins to mask the slow mode 
behavior. Positive noise on the measurement causes wit) 
to exceed u\ (<)V so that the logarithm becomes imaginary, 
resulting in gaps in the plot. 

The linearity of the response curve in region II is 
fundamental to the present technique, because it shows 
that there is only one spurious mode of concern with 
characteristic time in this range, as a second mode would 
cause extra terms w hich would create curvature. A slow - 
er mode would be of importance at larger times, which 
being in region III must be smaller than the noise level. 
A faster mode would appear at smaller times, that is, 
in region I. Such a mode would have appeared in testing 
of the point response function (Paden et al. 1997) il it 
were present. 

In order to plot Fig. 7, a value of vv. ls> is required. 
The initial estimate comes from inspection of Fig. 2. If 
the value used for vr as> is too large, the curve of Fig. 7 
will approach a constant value. If the value is too small, 
the curve w ill turn sharply down. By noting the behavior 
of ln|(vr. lsv ~ u’(/))/h’. 1s> - w nmi )| as a function of time, 
one can iterate on the value of vv as> to get its value quite 
accurately. Sufficient accuracy is achieved w hen the re- 
sponse is straight in region II and the initial part of 
region III. 

We now' consider the computation of the modal load- 
ing <\ If the input to the slow mode were a step function, 
then from Eq. (19); 


1 T W’as> ~~ l, ’mm 

where vr u is the measurement after the rapid response 
of the undistorted signal, that is. the primary mode of 
the detector, has damped out. The input to the slow mode 
is not a step response, but is the response of the detector 
to the source as the instrument scans onto the source, 
which constitutes region I. Because the primary mode 
of the detector is fast compared to the slow mode, that 
is, these two modes are well separated, region I is of 
short duration compared to region II. We select the mid- 
point of region I as the time at which a step input is 
applied. The straight line of region II is then extrapo- 
lated to this r 0 to give u () , so that the right-hand side 
of Eq. (21) and c can be computed. Fortunately, the 
slope of the line, -A, is small so that c is not sensitive 
to the estimate of / () . 

The numerical filter was applied to the PFM total 
channel calibration data and the results are shown in 
Fig. 2 as a dotted line using c = 0.026 and A = 9.45. 
The measurements approach the asymptotic value with 
a characteristic time of 9 ms, corresponding to the re- 
sponse time of the detector. The slow mode does not 
appear, indicating that the numerical filter does work 
and the parameters selected are suitable. 


5. Flight results of numerical filter for PFM 

The numerical filter was applied in the data process- 
ing for the CERES PFM instrument aboard the TRMM 
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Fig. 8. Average total channel radiances for I h as a function of 
sample number for PFM in flight. 

spacecraft and to the flight models I and 2 (FM-1 and 
-2) aboard the Terra spacecraft. The CERES instrument 
scans from the space look at 8° (horizontal is taken as 
zero) across the earth to another space look at 162° and 
then to 194° (i.e., 14° above horizontal), where it stops 
for a short time to view the Internal Calibration Module 
(Smith et al. 1998). The measurements from the PFM 
total channel for a scan on 5 August 1998 are shown 
as a function of sample number in Fig. 8. The radiances 
computed without the numerical filter are shown as a 
solid line and with the numerical filter by a dotted line. 
The same gains are used for both cases. The difference 
between the radiances with and without the numerical 
filter is also shown as a dashed line with the enlarged 
scale on the right axis. Each scan takes 6.6 s and a 
measurement is taken once every 1 0 ms so that at sample 
number 660 the scan cycle repeats. The scan begins with 
the space look at 8° scan angle to establish zero radiance. 
At approximately sample 60 the instrument begins to 
scan across the earth. At approximately sample 290 the 
instrument scans from earth to space (near scan angle 
162°) and at sample 310 scans to look at the Internal 
Calibration Module (ICM), where the total channel 
views the Internal Blackbody. At sample 350 the in- 
strument scans from the ICM to space (again near scan 
angle 162°) and at sample 380 has scanned onto the 
Earth for another scan, until sample number 610, where 
the instrument has returned to the starting position. The 
effects of the numerical filter are seen in regions w here 
there is a large change, that is, as the instrument scans 
onto the limb of the earth, as it scan onto the internal 
blackbody and as it scans to the space look. The rest of 
the time, the effect of the numerical filter is very small 
in the average. 

Figure 9 shows the approach of the total channel mea- 
surements to equilibrium while looking at the internal 
blackbody. Without the numerical filter (solid line) the 


Numerical Filter - OFF 
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FiCi. 9. Average total channel response for 1 h as it scans onto 
internal blackbody in High t . 


measurement does not reach steady state during the time 
available for the internal calibration and causes the mea- 
surement to be 0.5 W rrr 2 sr~' less than the asymptote 
at the end of the internal calibration period. The re- 
trieved undistorted measurements, that is, with the nu- 
merical filter (dotted line), quickly come within 0.1 W 
m 2 sr 1 of its asymptote, with the 9-ms response time 
of the detector apparent. If not accounted for, this dif- 
ference together with the difference at the space look 
would cause an error in the inflight calibration results 
for the gain of the channel. This error of gain would 
then affect all measurements. 

The output of the detector is given by the number of 
counts above the zero radiance level, w hich is converted 
to radiance units by multiplying by the gain of the chan- 
nel. The zero level drifts over a period of several minutes 
and is established for each scan by the space look. Figure 
10 shows the space look measurements in detail. As the 
instrument scans from viewing earth to space at scan 
angle 8°, the slow mode is obvious in the measurements 
without the numerical filter applied (solid line). The zero 
radiance value is taken to be the average of the mea- 
surements from 13 samples, typically sample number 
28 to 40. During the space look it is possible to get 
small negative values of radiance due to random errors 
below this average value. With the numerical filter ap- 
plied, the retrieved measurements very quickly drop as 
expected for the 9-ms first-order response of the detector 
and remain steady over the period of the space look. 
Any error in the space look causes an error in each 
measurement of earth radiances. No trace of the desir- 
able transient is discernible in the retrieved measure- 
ments; only residual noise at the 0.02 W m 2 sr -1 level 
remains. The space look is especially valuable for val- 
idating the numerical filter because it demonstrates that 
the transient mode is due to effects within the scanner 
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Sample Number 


Fi(i. 10. Average total channel counts for space look for 1 h as a 
function of sample number for PFM. 


alone and is not due to a mode which may be conjectured 
as created by a radiative coupling between the telescopes 
and the internal calibration sources when data from 
these devices are used. Also, the flatness of the retrieved 
measurements during the space look provides validation 
of the assumption of a single linear mode to describe 
the spurious transient. 

The radiances from the total channel of the PFM dur- 
ing the space looks just prior to and after the internal 
calibration are shown in Fig. 1 1. These space looks are 
not usually used in data processing, but are for vali- 
dating that this space look is compatible with those at 
scan angle 8°. Without the numerical filter, the radiance 
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Fid. 12. Shortwave channel response U> the SWICS for FM-1. 


measurements do not reach the steady state value. By 
applying the numerical filter, the retrieved radiances in- 
dicate the zero radiance level to be -0.06 W m 2 sr 1 , 
which would be a much smaller error than the 0.3 W 
m 2 sr 1 indicated by the raw r measurement. There are 
cases in which the space look is contaminated or lost 
in the data stream for various reasons. In such cases, 
the space look near scan angle 162° can be used after 
application of the numerical filter. 

The application of the numerical filter to the CERES 
Flight Models 1 and 2 aboard the Terra spacecraft is 
now described. Figure 12 shows the shortwave channel 
response to the Shortwave Internal Calibration Source 
(SWICS) for FM-1. There is a small overshoot of the 
response above the aymptotic value as the detector scans 
onto the SWICS and again as it scans away from the 
SWICS. due to extraneous reflections at the edge of the 
SWICS. Figure 13 shows ln|(u\ lsy — w(t))/(w. lsy - w,,)] 
for these measurements. The response overshoot results 
in the dip at sample number 192, after which the curve 
shows the linear behavior of the slow mode. The over- 
shoot was so quick that the slow mode does not sig- 
nificantly respond. From the slope of the curve, A is 
computed to be 2.9 s 1 and extrapolation of the line 
back to the initiation of the SWICS input gives c = 
0.016. 



Fk;. 1 F Average total channel response for I h for space look on 
internal calibration module side of scan for PFM. 


Fig. 13. ln[(vr K1 , - u(/))/(n lsv - uv,,,,, >] as a function of sample 
number for FM-1 shortwave channel as it views the SWICS. 
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Fig. 14. Total channel response to internal hlaekbody tor FM-2 
instrument with and without numerical tiller. 


Figure 14 shows the effect of applying the numerical 
filter to the total channel data from the FM-2 instrument 
during the use of the internal hlaekbody (IBB) for in- 
night calibration checks. The solid line is the response 
of the detector without accounting for the slow mode. 
The dashed line shows the response with the slow mode 
taken into account using the parameters determined 
from ground calibration data and the best fit parameters 
from in-orbit data are used for the dotted line. As the 
total channel scans to look at the IBB, it sees shortwave 
and longwave radiation scattered and emitted from the 
interior of the instrument, which causes the overshoot 
peak of 1 .6 W m 2 sr 1 . After the detector stops at the 
IBB position, the slow mode continues to vary through- 
out the duration of the IBB measurement. Between sam- 
ple numbers 320 and 343 the radiance value with the 
slow mode numerically filtered out varies less than 0.1 
W m 2 sr The parameters determined from ground 
calibration data and from flight data give results that are 
nearly indistinguishable. For this channel the effect of 
the numerical (liter is to change the radiance by 1.2 W 
m 2 sr ! at the beginning of the internal calibration 
period. 

Figure 15 shows the effects of the numerical filter on 
the space look radiance values for the total channel. The 
sample numbers start with zero at the center of the space 
look, after which the detectors scan the earth, take a 
short space look on the other side, and then look at the 
calibration devices, scan back across the earth to the 
initial space look and repeat (Smith et al. 1998). The 
abscissa begins with the scanning of the detectors to 
space and continues across the zero sample number po- 
sition at which point a new' data frame begins. The ra- 
diance for each position is the mean for that sample 
number for scans taken over a 1/2-h period. With the 
numerical filter, the radiance level is 0.04 W m 2 sr 1 
by sample number 615. but without the numerical filter 
the radiance level does not drop to 0.04 until sample 
number 652. Application of the numerical filter consid- 
erably increases the length of time available for the 
determination of the zero radiance level. This is espe- 
cially important during times when artifacts or data dis- 


Wlthout Filter 

With Filter, C = 0.025, X = 9.20 (Flight Data) 



Fig. IS. Tola! channel response during space look tor FM-2 
instrument with and without numerical filter. 


ruptions occur in the space look. One such case is when 
the moon appears in the space look. 

The shortwave channel response for the space look 
is shown by Fig. 16 for FM-2. The values of the pa- 
rameters computed from ground calibration (c = 0.023, 
A = 8.2) result in an undershoot of the radiance but 
parameters (c = 0.013, A = 8.6) computed using flight 
data and a best mean square fit procedure give radiance 
values which are within 0. 1 W m 2 sr 1 of 0 for sample 
number 612. The flight parameters will be incorporated 
into the data processing for the FM-2 shortwave chan- 
nel. Additional information about the FM-1 and -2 slow 
mode characterizations, including the description of the 
best mean-square fit method is given by Smith et al. 
( 2000 ). 

6. Discussion 

As the instrument scans from the space look to the 
earth, the effect of the slow mode is to cause the mea- 


Wlthout Flits r 

With Fitter,- C - 0.013, X « 8.60 (Flight Data) 
With Filter. C = 0.023, X « 8.21 (Ground Data) 



Sample Number 

Fig. 16. Shortwave channel response for space look for FM-2 
instrument with and without numerical filter. 
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surements to be low until the instrument has scanned 
through 40° of nadir scan angle, by which time the effect 
of the slow mode has decreased. The numerical hi ter 
eliminates this effect. When the instrument scans from 
one scene to another with greatly differing radiances, 
for example, from tropical land to the top of a deep 
convective cloud, the instrument indicates an OLR ra- 
diance which is high initially by approximately 2% of 
the change, due to the slow mode. As the instrument 
scans from the cold cloud top to warm land, the lag 
causes the instrument to indicate a radiance which is 
low initially. These differences average out and do not 
appear in averages over large areas that are well away 
from the limb. For this reason, the major effect of the 
numerical filter is in establishing the correct zero ra- 
diance value during the space look and the correct re- 
sults from the internal calibration data. These values 
affect all measurements. 

The fact that the spurious signal is small and its time 
response is well separated from the rest of the system 
response modes simplified the implementation of the 
numerical filter, but is not a requirement in the design 
of the numerical filter. Likewise, if additional modes 
were present in the instrument output, they would com- 
plicate the development a numerical filter, but would 
not prevent it. 

Further resolution enhancement of the image is not 
attempted with this numerical filter, as the CERES in- 
strument w r as designed to minimize the sum of blur and 
alias errors near nadir. It may be possible to design a 
numerical filter that would enhance resolution in the 
along-track direction away from nadir, where the over- 
lapping of pixels increases greatly, but that is beyond 
the intent of this paper. 

7. Conclusions 

During calibration testing of the CERES proto-flight 
model, a spurious slow response was found with a char- 
acteristic time of 0.3-0.4 s, an order of magnitude slow- 
er than the primary time response mode. A numerical 
filter is developed here for removal of this spurious 
response. This method does not depend on the details 
of the underlying physics of the transient, but only re- 
quires that the response be linear. It is demonstrated that 
the slow response is due to a single mode. Ground cal- 
ibration data are used to compute the modal response 
time and amplitude for the spurious mode. Application 
to flight data for the proto-flight model and for flight 
models 1 and 2 aboard the Terra spacecraft shows that 
the numerical fi Iter eliminates the slow mode from the 
measurements; only residual noise from other effects 
remain. 

The effects of the slow : mode are primarily in the in- 
flight calibration of the instrument, whereby the slow 
mode prevents the instrument from attaining a steady 
state during the time available for the space look, which 
is required to establish the zero-radiance counts and for 
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the standard in-flight calibration sequence. The numer- 
ical Filter does not affect averages over large areas aw'ay 
from the limb of the earth onto which the instrument 
scans from space. 
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APPENDIX 

Relation of Spurious Transients to Undistorted 
Measurements 

The spurious transients in the measurements are due 
to slow' modes of heat translet* within the thermistor 
bolometer assembly. The relation ol these transients to 
the undistorted measurement is based on the time-de- 
pendent temperature distribution within the assembly, 
which consists of the thermistor bolometer, the paint 
layer which absorbs the radiation to be detected and 
which is bonded to the thermistor bolometer, and the 
active and compensator disks. The equation for the tem- 
perature distribution within the detector assembly is 

c—T= V- (kVT) + q. ( A 1 ) 

(H 

where q is a heat source in the detector. The heat source 
of primary interest is the radiation that is absorbed in 
the layer of paint that is bonded to the thermistor bo- 
lometer. The method of separation of variables is used, 
whereby it is assumed that T{ r, t) — X(r )</>(/), whence 

-A f ,cX n = V • (kVX„) + q„ and 

— + A >„ = (A2) 

dt 

where q n is the loading of the heating in the n - th mode 
and /„(/) is its time variation. The temperature at a point 
is then 

T( r. /) = 2 X,(r)< P„U). 

where C denotes the set of modes which contribute to 
the temperature. One can include electrical feedback 
modes in the analysis provided they are linear. Their 
addition simply adds more equations to the set of Eq. 
(A2). The measurement of the thermistor bolometer is 
the integral of the temperature over the region between 
the electrodes and thus can be written as 

vi' (f) - 2 /?„<£„(»• (A3) 

>u ( 

The complete set of modes C includes the primary re- 
sponse time of the detector, a number of very fast tem- 
perature modes to which the thermistor bolometer has 
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very little response, and also four inodes for the re- 
sponse of the Bessel filter. These modes determine the 
undistorted signal u(t). The set C also includes a set S 
of slower modes which are spurious and are to be taken 
into account by the numerical filter. Equation (A3) thus 
becomes 


w(t) - i4(t) + u(0, (A4) 

where v(t) is the spurious part of the signal: 

V(l) = 2 b„4> n (t). (A5) 

/ig .S' 


We now make the approximation, which is funda- 
mental to the application of this method: For the slower 
modes, the radiative forcing f„(t) is proportional to the 
undistorted signal u(t). This approximation relies on the 
u( 0 response being fast relative to each of the slow 
modes. 

The operator Y\ neS (d/dt + A„) is applied to each side 
of Eq. (A5) and the approximation is used to replace 
/„(/) with u(t) times a constant for each n , resulting in 



S(n) 


s c, n (^ + a /.) w > ( A6 ) 


and S(n) denotes the set S of spurious modes with the 
nth mode removed. The c n is the b„ times the propor- 
tionality of response of the nth mode to «(/). The left- 
hand side of Eq. (A6) is a polynomial in the time de- 
rivative operator dfdt, of degree equal to the number of 
spurious modes which must be taken into account; the 
right hand side is a polynomial of one degree less. Equa- 
tion (A6) can be written as 


8 



v(0 



(A7) 


This equation relates the spurious modes to the undis- 
torted measurements. The solution of Eq. (A7) is 


u(0 = 2 


a n + I u(t)c jAt dr 


However, u(t) is not observable. Although this equation 
was derived on the basis of time-dependent heat transfer, 
the same form of the equation would be found, for ex- 
ample, if the spurious transients were due to time-de- 
pendent electronic responses, provided only that the 
mechanism is linear. 
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